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Abstract 

The radiative instability of a beam moving in a photonic crystal of finite dimensions 
is studied. The dispersion equation is obtained. The law r ~ pVC'+a) is shown to 
be valid and caused by the mixing of the electromagnetic field modes in the finite 
volume due to the periodic disturbance from the photonic crystal. 



Numerous works analyze the generation of induced radiation by electron beams 
moving in spatially periodic media [T|2|l3|l^|f5] . Thoroughly studied are the pro- 
cesses of wave generation in a one- dimensional case, when the electron beam 
moves along the axis of a corrugated waveguide (traveling- wave tube (TWT), 
backward- wave tube (BWT)) or along the axis of the undulator (volume free 
electron lasers, ubitrons). It was found out that under the action of radia- 
tion, the electron beam with a uniform density distribution becomes spatially 
modulated, i.e., the radiative instability of a beam emerges. The increment of 
the electron beam instability is a most important quantity characterizing the 
generation capability of the beam. The analysis of all mechanisms of induced 
radiation generation by relativistic beams in the Compton regime showed that 
the increment of beam instability F for a cold beam (i.e., a beam where all 
the electrons have the same longitudinal velocity u) follows the law Y ~ p 1 ^ 3 , 
where p - is the beam density [2|3] . 



In 1984 it was shown [6] that induced X-ray radiation produced by the elec- 
tron beam passing through the crystal under the conditions providing the 
coincidence of roots of the dispersion equation, which describes the relation 
between the wave vector k and the photon frequency 00 in the crystal, leads 
to the appearance of a new law for the increment of the beam instability: 
T ~ p 1 /( s + 3 ) ; where s is the number of additional waves appearing due to 
diffraction of emitted X-ray quanta in the crystal. This law of instability leads 
to a significant reduction of the laser generation threshold for X-ray radiation 
in the crystal (according to this law, the generation threshold for a LiH crystal 
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can be achieved when the beam current density is 10 s A/cm 2 in contrast to 
10 13 A/cm 2 as required according to the conventional law p 1 / 3 . 

It was later shown in [7|8] that the above described law is valuid for all wave- 
length ranges of induced radiation produced by electrons in spatially periodic 
media (diffraction gratings) and for various types of nonlinear interaction of 
waves in both natural and artificial (or frequently called "photonic") crystals. 
The results obtained in these works enable design and development of a new 
type of Free Electron Lasers called the Volume free Electron Lasers (VFEL) 



Theoretical description of the generation processes in a photonic crystal placed 
inside the resonator was given in [9|12|I13] . The first and most important step in 
describing the generation process in VFELs (FELs and so on) is the analysis of 
the problem of the electron beam instability in the resonator. The theoretical 
study of the instability of electron beams moving in natural and artificial 
(photonic) crystals was carried out for the ideal case of an infinite medium 
(see the review in [9] and Pll0fllfl2|[T3] ). The question arising in this regard 
is how the finite dimensions of the photonic crystal placed inside the resonator 
affect the law of electron beam instability. It is known, for example, that the 
discrete structure of the modes in waveguides and resonators is crucial for 
effective generation in the microwave range [2"f5flfo'] . 

In present paper the radiative instability of a beam moving in a photonic 
crystal is studied. The dispersion equation describing instability in this case 
is obtained. It is shown that the law T also valid and caused by 

the mixing of the electromagnetic field modes in the finite volume due to the 
periodic disturbance from the photonic crystal. 

The system of equations describing generation of induced radiation in photonic 
(and natural) crystals can be obtained from Maxwell equations: 



-* 1 3D 4tt- - ldH 

r0tH =clH + -^ r0tE = -clK> (1) 

-> dp -> 

divD = 4np, — + divj = 0, 



— * — * 

where E and H are the strength of the electric and the magnetic field, respec- 
tively; j p are the current and charge densities; Di (r, t') = J Su (r, t — t') E\ (r, t') dt' 
or Di(f,u) = Eu (r,u) E[ (r,u), where indices i,l = 1,2,3 correspond to 
x, y, z; En (r, u) is the dielectric permittivity tensor of the photonic crystal. 
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The current and charge densities are defined as: 

J(r,t) = e^vjft- r a {t)), pft,t) = e^<5(f- f a (t)), 



ft) 



where e is the electron charge, v a is the velocity of the electron with number 
a in the electron beam, 



-£ = \E(f a (t),t) + - 

at m^ a I c 



(t) x H (r a (t) ,t)]-^ (v a (t) E (f a (t) ,t))\, (3) 



where j a = (l — ^\ 2 is the Lorentz factor, E ft a (t) , t) (H ft a (t) ,t)) is the 
electric (magnetic) field strength at point f a , where the electron with number 
a is located. Note that equation ([3]) can be written as follows 



dp a 

dt 



m- 



djaVa 

dt 



E (f a (t) ,t) + - v a (t)xH (f a (t) , t) 
c 1 



(4) 



where p a is the particle momentum. 
From equations (TS]) one can obtain 

A - - / - 1 d 2 D An dj 



The dielectric permittivity tensor can be presented in the form e (r) = 1+x ft), 
where % ft) is the susceptibility. 

For x « 1, equation (|3j) can be rewritten as 

^ ft V-^fz ft t-^E ft, dt' = 4tt (^^J^ + Vp ft, *)) .(6) 



In the general case, the susceptibility of the photonic crystal reads xft = 
J2 Xce.ii ft ' — ity, where Xceii ft — ri) is the susceptibility of the crystal unit cell. 

i 

The susceptibility of an infinite perfect crystal x ft can be expanded into the 
Fourier series as follows: x ft = J2Xr^ Tr , where f is the reciprocal lattice 

T 

vector of the crystal. 

To be more specific, let us consider in details a practically important case 
when a photonic crystal is placed inside a smooth waveguide of rectangular 
cross-section. 
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The eigenfunctions and the eigenvalues of such a waveguide are well-studied 

USEE]. 



Suppose the z-axis to be directed along the waveguide axis. Make the Fourier 
transform of ((Sj) over time and longitudinal coordinate z . Expanding thus 
obtained equation for the field E (f±, k z ,u) over a full set of vector eigenfunc- 
tions of a rectangular waveguide Y^ n (r±, k z ) (where m,n = 1,2, 3...., while A 
describes the type of the wave [T7] , one can obtain for the field E the equality 

E (f ± , k„ uj) = Y: C x mn (K w) YL (r±,h) . (7) 



As a result, the following equation can be written 



(kg + K^nx) Cmn (kz j w ) 

-7S E / ^ fc) X (r) Y* v (r ± ,k' z ) e -(^)* d 2 r ± C* n , (k' z , u) dk' z dz = (8) 



m'n'X' 



Aniui 



I Y^n (r±, k g ) {j(r±,z,uj) + (Vj(r ± , z, e lk * z d 2 r ± dz 



where ft^ nA — + ^ynA- 

The beam current and density appearing on the right-hand side of flS]) are 
complicated functions of the field E. To study the problem of the system 
instability, it is sufficient to consider the system in the approximation linear 
over perturbation, i.e., one can expand the expressions for j and p over the 
field amplitude E and abridge oneself with the linear approximation. 

As a result, a closed system of equations comes out. For further consideration, 
one should obtain the expressions for the corrections Sj and Sp due to beam 
perturbation by the field. Considering the Fourier transforms of the current 
density and the beam charge j (k, ujj and p (k, cSj , one can obtain from (j2J) 
that 

— * / — * \ { , -. s. kSv a (u -ku a )) 
Sj (k, u) = e e" \ Sv a [u - ku a ) + u a '- \ , (9) 

a=l [ k> KU a J 



where f ao is the original coordinate of the electron, u a is the unperturbed 
velocity of the electron. 

For simplicity, let us consider a cold beam, for which u a ~ u, where u is the 
mean velocity of the beam. The general case of a hot beam is obtained by 
averaging Sj (k, u) over the velocity u a distribution in the beam. 
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According to (jSJ), the velocity Sv a is determined by the field E(r a ,u) taken 
at the electron location point f a . The Fourier transform of the field E (f a ,u) 
has a form 

E(f a ,u) = ^3 J E(k',uj) e^d'k'. 



As a result, the formula for 5j (k,u) includes the sum J2 e ^ k ) ra over the 

particle distribution in the beam. Suppose that the electrons in an unperturbed 
beam are uniformly distributed over the area occupied by the beam. Therefore 



£ e -i(k-k') r a = ^3 Pq5 ^_ # 



where po is the beam density (the number of electrons per 1 cm 3 ). 

As a result, the following expression for Sj (k,ojj can be obtained [T8fT9] : 

\ iue 2 p (k 2 - ^) -> \ 

Sj (k, u) = ^ 2 c > uE (k, u>) . (10) 

iui — ku) m^juj 

Using the continuity equation, one immediately obtains the expression for 
p (k,ujj. Expression (fTUj) . the inverse Fourier transform of E (k,uj, and the 
expansion (j7j) enable writing the system of equations (jSj) as follows: 

{k 2 z + H 2 nn \) — ^ C^n (^2!^) — 

-Si E / (f±,k,) X (r) i^v £u e-<*--^ rfVC*' v (tf„ w) dtf/fe = (n) 

J^i (£]_,&*) d 2 £af} C* n , 



m'n'A' 



7C 4 (w-fc mn «)^ 



where k z ) = J e"^^^ (r ± , rf 2 r ± . 

Note that within the limit where the transverse dimensions of a photonic 
crystal tend to infinity, the expression between the braces takes the form 
(eu) 2 , where e is the unit polarization vector of the wave emitted by the beam 
[T81IT9]. 

Now let us consider the integrals on the left-hand side of equation (jlip . 
Note that according to [T5lll6|ll7] . the eigenfunctions Y^ n (f±,k z ) of a rect- 
angular waveguide include the combinations of sines and cosines of the form 
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sin^x, cos^x (sin^y , cos^y) , i.e., in fact, the combinations e % ™ x ,e l ™ y . 
Hence, the left-hand side of the equation includes the integrals of the type 

I=j e ~ l ^ x Xceii -Xi,y- y u z - z { ) e l ^ x dx. 



The substitution of variables x — X{ = rj gives the sums of the form 

i 

where Xi = d x fi, where d x is the period of the photonic crystal along the 
x-axis, fx = 1, 2, ...N x , where N x is the number of cells along the x-axis. 

The above-mentioned sum 

■jr(m—m')d x N x 

$ e -il{m-m')xi _ e i£(m-m'KN x -l)ds SlU 2a tyft 

x /-^ • ir(m—m')d x ' \ ' 



If m — mf — 0, then S x = N, 



X ■ 



Let us now discuss what this sum is equal to when m — m! = 1. In the 
numerator d x N x = a, hence, the nominator is equal to 1 (sin^ = 1), while in 
the denominator sin 1 ^- rs As a result, the relation Sx( f n ~ m -*) — 1 ~ 0.6. 

2a 2N X ' S x (m— m — 0) 7r 

With growing difference m — m', the contribution to the sum of the next terms 
diminishes until the following equality is fulfilled 

= nP, (13) 



where P = ±1, ±2... In this case the sum S x = N x . 

The similar reasoning is valid for summation along the y-axis. 

It follows from the aforesaid that if the equalities like (fT2]) . (fl3"|) are fulfilled, 
that is, the equalities k xm — k' xm , = t x are fulfilled (i.e., k' xm , = k xm — r x ), where 
t x = ^f-F is the x-component of the reciprocal lattice vector of the photonic 
crystal, F — 0, ±1, ±2... and k yn — k' yn , = r y (i.., k' yn , = k yn — r y ), where r y = 
^■F' is the ^-component of the reciprocal lattice vector of the photonic crystal, 
F' = 0,±1,±2...), then the major contribution to the sums comes from the 
amplitudes C*' v (k' z , u) = C x ' (k ±mn -T±,k z -r x ,u) = C x> (k mn - f, u). 

In describing the system we shall further consider only those modes that satisfy 
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the equalities of the type ( !T2|) . (jl~3j) . As stated above, the contribution of other 
modes is suppressed. 

As a result, one can rewrite the system of equations (fTTj) as 



(k mn — ^rj C X (k mn , U?j — J2 Xmn (t)C X (kmn T, LO 

JuY x n (k±,k z ) d 2 k ± \ 2 \c x (k mn ,u) , 



2tt 



(14) 



i.e., 



1 "I - Xmn (0) 



^Ex^(r)C A ' (k-r,u)=0 

X't V 7 



2- 



(15) 



where X™ (r) = £ / F„„ (r ± , k z ) \ f±, T z )Y*, n , (f ± , k z - r z ) d 2 r ±1 \ fx, r z ) = 
I Xceii (x — Xi, y — i/i, () e~ lTz< *d(, m! and ml are found by the conditions like 

f lT3|) . ujl is the Langmuir frequency, oo 2 = Ane m p0 ■ 

This system of equations coincides in form with that describing the instability 
of a beam passing through an infinite crystal [l"8l[T9] . The difference between 
them is that the coefficients appearing in these equations are defined differ- 
ently and that in the case of an infinite crystal, the wave vectors k mn have a 
continuous spectrum of eigenvalues rather than a discrete one. 

These equations enable one to define the dependence k (co), thus defining the 
expressions for the waves propagating in the crystal. By matching the incident 
wave packet and the set of waves propagating inside the photonic crystal using 
the boundary conditions, one can obtain the explicit expression describing the 
solution of the considered equations. 

The result obtained is formally analogous to that given in [20] . 

According to (|15|) . the expression between the square brackets acts as the 
dielectric permittivity e of the crystal under the conditions when diffraction 
can be neglected: 



e = n 2 = l + Xm \ (0) 



U L (krnn c2 



'7c 2 (co - k mn uj 



2 1 2tt 



a 



n is the refractive index. 
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As is seen, in this case the contribution to the refractive index comes not only 
from the scattering of waves by the unit cell of the crystal lattice, but also 
from the scattering of waves by the beam electrons (the term proportional to 
the photonic crystal penetrated by a beam of electrons is a medium that 
can be described by a ceratin refractive index n (or the dielectric permittivity 
eo)- 

According to (IT51) . the beam contribution increases when u — > ku. 

Since this system of equations is homogeneous, its solvability condition is the 
vanishing of the system determinant. 

In the beginning, let us assume that the diffraction conditions are not fulfilled. 
Then the amplitudes of diffracted waves are small. In this case the sum over 
r can be dropped, and the conditions for the occurrence of the wave in the 
system is obtained by the requirement that the expression between the square 
brackets equal zero. 

This expression can be written in the form (the velocity {t||oz) 




where n is the refractive index of the photonic crystal in the absence of the 
beam e = n\ = 1 + x„„ (0), 

i.e., 




Since the nonlinearity is insignificant, let us consider as the zero approximation 
the spectrum of the waves of equation ffl6|) with zero right-hand side. 

Let us concern with the case when u — k z u — > (i.e., the Cherenkov radiation 

condition can be fulfilled) and (k 2 z — (^^o — K mn)) ~~ * 0^ i- e 5 the electromag- 
netic wave can propagate in a photonic crystal without the beam. With zero 
right-hand side the equation reads 



CO 



_ 2 _ 2 

2 n K mn 



(00 



k z u) = 



(17) 
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As a consequence, in this case the roots of the equation are 
2 _ ^mn_ J-' - _jk K 



c 



Since fc2 2 = - > in view of the Cherenkov condition, we are concerned with 
the propagation of the wave with k\ z > in the photonic crystal. In this case 
in the equation for k z , one can take (k z — k\ z ) (k z + k\ z ) ~ 2k\ z (k z — k\ z ) and 
rewrite equation ( fl6l) as follows: 



(k z - hz) (k z - k 2z y 



uj 2 l u 2 (n 2 



2k lz u 2 ~/c 4 { 2ir 



uY^ n (k±, k z ^ d 2 k A 



(19) 



i.e., 



(kg - k u ) (k z - k 2z f = -A 



(20) 



where A is real and A > (as for the occurrence of the Cherenkov effect, it 
is necessary that Uq > 1). We have obtained the cubic equation for k z . Let us 
consider the case when the roots k\ z and k 2z coincide k\ z = k 2z . It is possible 
when the particle velocity satisfies the condition 



u 



(21) 



Tin - 



Introduction of £ = k — k\ z gives for k\ z = k 2z 

e = -a. 



(22) 



The solution of equation fl22|) gives three roots £i = —\/A, £2,3 = \ (l ± v^- 

As a consequence, the state corresponding to the root £ 2 = \ fl + iV$J \f~A 
grows with growing z, which indicates the presence of instability in a beam 
[2Tj . In this case Imk z = Im^ 2 ~ &~p. 

Note here that the photonic crystal built from metallic threads has the re- 
fractive index tiq < 1 for a wave with the electric polarizability parallel to the 
threads, i.e., in this case the Cherenkov instability of the beam does not exist 
[9] (but if the electric vector of the wave is orthogonal to the metallic threads, 
the refractive index is no > 1 , so for such a wave the Cherenkov instability 
exists [22]). 
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It should be pointed out, however, that, unlike an infinite photonic crystal, 
the field in the crystal placed into the waveguide has a mode character, and so 
the presence of in the denominator of equation ( l2"Tj) results in reduction 
of the radicand in (12 ip to the magnitude smaller than unity even when > 1. 
Hence, u > c, which is impossible. Consequently, the radiative instability of 

2 2 

the above type in the waveguide can arise under the condition Uq — Krn - n 2 > 1 
rather than n 2 , > 1. 



Suppose now that in the photonic crystal the conditions can be realized under 
which the wave amplitude C mn (k mn + rj is comparable with the amplitude 

Cmn (kmn)- By analogy with the standard diffraction theory for an infinite 
crystal |23II24] . in the case under consideration, when \ << 1> h is sufficient 
that only the equations for these amplitudes remain in f[T5]) . 



To be specific, let us further consider a photonic crystal formed by parallel 
threads. Also assume that they are parallel to the waveguide boundary (y, z). 

Analysis of diffraction of a A-type wave with the electric vector in the plane 
(y, z) (a TM-wave) gives 



k 2 -— e 



C x (k mn ,uj)-^x^ n (-r)C x (k r , 







(23) 



~2 



C x (k mn +t,u>)- —xtl if) C x (k mn , u) 



Since the term containing (cu — (k + r) uj is small when (u — kuj vanishes, 
in the second equation it is dropped. 



The dispersion equation defining the relation between k z and u is obtained by 
equating to zero the determinant of the system (1231) and has a form: 



(u — k z u] 



(kmn l? £ o) ((kmn + 

/ uY x n (k ± , k g ) d 2 k ± \ Z I (k 2 mn c 2 - u 2 ) ( (k mn + r)' 



l)± 

4 1 2vr 



(24) 



Because the right-hand side of the equation is small, one can again seek the 
solution near the points where the right-hand side is zero that corresponds 
the condition of occurrence of the Cherenkov radiation and excitation of the 
wave which can propagate in the waveguide: 

(k 2 - (^6 - K 2 mn )) ({k z + T) 2 - (^E - (K mn + T ± ) 2 )) - ^XrX-r = 
2 



- s) 
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The roots of equations are sought near the conditions k 2 mn (k mn + t) , 



2 n 2 



k z = k z0 + £, k 2 z = k 2 z0 + 2k z0 C + e, k 2 z0 = ^e - K 2 mn , kz0 = ^s -^-(2Q) 

(K + T z ) 2 = [(k 0z + T z ) + if = (k 0z + T z f + 2 (k 0z + T 2 )£ + ^ 2 

Hence, 

{k 0z + r z f + {K rnn + f ± ) 2 + 2 (A; 02 + r z ) + 2 (A; 02 + r 2 ) £ + £ 2 = ^ 

(k mn + f) 2 + 2 (fc 02 + T Z ) e + ^ = fcLn + 2 ^ m nT + T 2 + 2 (Afo + T 2 ) £ + 

And one can get 

2k 0z £ (2 (A: 02 + T z ) £ + (2fc 0m nT + r 2 )) - — X rX-r = 



4k 0z (k 0z + t z ) i 2 + 2A: 02 (2k 0mn T + r 2 ) f - — XrX-r = (28) 



c 4 



^ + ^ L£ _ Z = 

(k 0z + T 2 ) c 4 4/c 02 + Tz) 



(2k r + 

6,2 = - V w , , ( ± 



T 2 



2 

2k T + T 2 \ UJ A XtX-t 



4 (k 0z + t z ) \J ^4 (& 02 + r 2 ) y c 4 4/c 02 (A*,* + r 2 ) 

If (k 0z + r 2 ) = — |fco z + r z \, the root can cross the zero point. At the same 
time, the second equation should hold. 

oo — k z u = oj — k 0z u — £u = 0. 
Consequently, 



u- k 0z u _u _u u K 2 mn c 2 

K — — ft o 2 — \\ £o y~ ■ 

u u u c V w z 



If 5 < 1, then i = a (l - /^e - > 0, £ = £ - A; 

Let the roots £i and £ 2 coincide (6 = £ 2 )- This is possible at point 

2k T + T 2 CO 2 ^X7X~r 



^0z 



4(k 0z + r z ) c 2 ^Akozlkoz + r^ 
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here k 0z + t z < 0. 

The roots coincide when the following equality is fulfilled 

U C ^k 0z \k 0z + T z \ 



i.e., 

- = T t / and hz = -\ £o 5— 

u c ^o.l^ + r,! c V 



Let e < 1) then - > k 0z (since u < c), the situation for the solution - = 
koz — 1 ^SS EI — gets complicated and the Vavilov-Cherenkov condition is 

c y/ 4k 0z \koz+T z \ 

not fulfilled. 



Now let us consider the solution - = k Qz + " vx^x-r — _ At r 2 < the 

" c y/4k 0z \k 0z +r z \ 

difference k$ z + r z can be reduced so that the sum on the right would appear 
to become equal to -, and so one could obtain 4 coinciding roots. 

Interestingly enough, for backward diffraction, which is a typical case of fre- 
quently used one-dimensional generators with a corrugated metal waveguide 
(the traveling- wave tube, the backward- wave tube), such a coincidence of roots 
is impossible. 

Indeed, let the roots £1 and £2 coincide. In this case for the backward Bragg 
diffraction \t z \ ~ 2ko z , t z < 0. Then by substituting the expressions for 

k 0z = ^ \J Eq — ^ m j2 and Eq = = 1 + Xmn (0) an d retaining the first-order 

infinitesimal terms, the relation - m kn z + ^J^- can be reduced to the form 

_ + ±M\ < h i.e., the equality does not hold and 

the four-fold degeneracy is impossible. Only ordinary three-fold degeneration 
is possible. 

However, if Eq > 1 and is appreciably large, then in a one-dimensional case, 
the four-fold degeneracy of roots is also possible in a finite photonic crystal"*"! 

Thus, the left-hand side of equation (|24|) has four roots £1, £2, and a double 



The authors are grateful to K. Batrakov, who drew our attention to the fact 
that for an infinite crystal with Eq > 1, the intersection of roots is possible in a 
one-dimensional case. 
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degenerated root £3. Hence, equation fl2l]) can be written as follows: 



If the roots coincide (£1 = £2 = £3), one obtains (£ — £1) = £>,, i.e., £ — £1 = 
^B. 

The fourth root of B has imaginary solutions depending on the beam density 
as Imk z ~ p l Q A (the parameter B ~ w|, i.e., S ~ p , see the right-hand side 
of (j24p ). This increment is larger than the one we obtained for the case of the 
three-fold degeneracy. 

The analysis shows that with increasing number of diffracted waves, the law 

established in [6|T7||8] is valid: the instability increment appears to be propor- 

1 

tional to p s + 3 , where s is the number of waves emerging through diffraction. 
As a result, the abrupt decrease in the threshold generation current also re- 
mains in this case (the threshold generation current j t h ~ ttt^Tu — ft^j where 
L is the length of the interaction area). 

It is interesting that according to [12], for a photonic crystal made from metal- 
lic threads, the coefficients x ( T )> defining the threshold current and the growth 
of the beam instability, are practically independent on r up to the terahertz 
range of frequencies because the diameter of the thread can easily be made 
much smaller than the wavelength. That is why photonic crystals with the pe- 
riod of about 1 mm can be used for lasing in terahertz range at high harmonics 
(for example, photonic crystal with 3 mm period provides the frequency of the 
tenth harmonic of about 1 terahertz (A=300 micron). 

The analysis of laser generation in VFEL with a photonic crystal when the 
beam moves in an undulator (electromagnetic wave) located in a finite crystal, 
made similarly to the above analysis, shows that in this case the dispersion 
equation and the law of instability also have the same form as in the case 
of an infinite crystal. The procedure for going from the dispersion equations 
describing instability in the infinite case ( TT5|) is similar to that discussed earlier 
in this paper. It consists in replacing the continuous k by the quantified values 
of k mn and redefining the coefficients appearing in equations like (fT5]) . 

It is important to emphasize the general character of the rules found in this 
paper for obtaining dispersion equations that describe the radiative instability 
of the electron beam in a finite photonic crystal. In particular, they are valid 
for describing the processes of instability of an electromagnetic wave in finite 
nonlinear photonic crystals. 
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